Observation of Hysteretic Transport Due to Dynamic Nuclear Spin Polarization 

in a GaAs Lateral Double Quantum Dot 
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We report a new transport feature in a GaAs lateral double quantum dot that emerges for 
magnetic-field sweeps and shows hysteresis due to dynamic nuclear spin polarization (DNP). This 
DNP signal appears in the Coulomb blockade regime by virtue of the finite inter-dot tunnel coupling 
and originates from the crossing between ground levels of the spin triplet and singlet extensively 
used for nuclear spin manipulations in pulsed gate experiments. The magnetic-field dependence 
of the current level is suggestive of unbalanced DNP between the two dots, which opens up the 
possibility of controlling electron and nuclear spin states via DC transport. 

PACS numbers: 73.63.Kv, 73.21. La, 76.70.Fz, 72.25.Rb 
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Electron spin dynamics in a semiconductor quantum 
dot (QD) is affected by the nuclear spin (NS) ensemble 
of the host material through the effective magnetic field 
(Overhauser field) produced via the hyperfine (HF) in- 
teraction [iHEI. This is important in view of quantum 
information processing using physical qubits based on 
electron spins in QDs |l|-|3|. In particular, when a qubit 
is encoded in the Sz = subspace of singlet and triplet 
states in a double QD (DQD), even a small difference in 
the Overhauser field ABn = Bn,l — Bn,r between the 
left and right dots can drastically affect the time evolu- 
tion of the qubit, with its coherence time limited by the 
fluctuation in ABn- On the other hand, electron-nuclear 
spin coupled dynamics that takes place around a spin 
triplet-singlet (T-S) crossing has been exploited as a re- 
source for preparing both electron and nuclear spin states 
in a controlled manner Using various pulsed gate 

sequences, |Ai3^| has been enhanced or its fluctuation 
squeezed, which respectively allow for rapid and univer- 
sal control [7[ or a longer coherence time of the 
qubit. Such ambivalent roles of NSs have motivated nu- 
merous studies aimed at understanding and controlling 
NS dynamics in DQDs 

DC transport in the spin blockade (SB) regime [isj , 
where the charge transport is governed by transitions be- 
tween different spin states, is one of the powerful means 
for investigating the interplay between electron and nu- 
clear spins in DQDs. Experiments have revealed cur- 
rent fluctuations arising from random dynamics of NSs 
14, [T^ and/or hysteresis due to dynamic nuclear po- 
larization (DNP) However, since DC trans- 
port relies on sequential tunnel processes involving dif- 
ferent charge states, such as (m,n) — (m -I- 1, n — 1) — 
{m,n — 1) — > {m,n) [where {m,n) denotes the occupa- 
tion of the left and right dots], previous studies have 
been mainly focused on examining cases where all rel- 
evant electronic states are within the transport window. 
Consequently, electron-nuclear spin dynamics at the T- 
S crossing, which has been extensively used in pulsed 



gate experiments [6|-l8|, has been poorly explored in DC 
regime. 

In this Letter, we report a new transport feature in the 
Coulomb blockade (CB) regime of a few-electron GaAs 
lateral DQD. We observe a sudden increase in the dot 
current I^^t in the CB regime near the (1, l)-(2, 0) charge 
boundary. The feature shows hysteresis indicative of dot 
state locking due to DNP associated with the crossing be- 
tween the ground levels of the triplet and singlet. As the 
inter-dot tunnel coupling is reduced, the feature weakens 
and eventually disappears, which indicates that the hy- 
bridization between the (1, 1) and (2,0) singlets is essen- 
tial for the lifting of the CB. Magnetic- field dependence 
of /dot suggests that the DNP generates strong imbal- 
ance of NS polarization between the two dots, which our 
model shows to be possible via positive feedback mech- 
anism. Our results provide new insights into electron- 
nuclear spin dynamics in a DQD and open up the possi- 
bility of controlling electron and nuclear spin states via 
DC transport. 

As shown in Fig. 1(a), the DQD used in this study was 
defined with Ti/Au gates in a GaAs/Alo.sGao.yAs het- 
erostructure containing a two-dimensional electron gas 
(density 2.2 x lO^^m-^ and mobility 200mVVs) 80 nm 
below the wafer surface. We applied a fixed source-drain 
bias VsD = — 800/LtV and measured /dot as a function of 
gate voltages Vrl and Vpr and in-plane magnetic field B\\ 
applied in the direction shown by the solid arrow in the 
figure. The energy levels in the two dots were controlled 
with VpL and Vpr, while the inter-dot tunnel coupling 
t was tuned with Vc- The actual values of (rn,n) were 
determined using a side-coupled quantum point contact 
charge sensor. All measurements were carried out with 
the sample mounted in a dilution refrigerator with a base 
temperature of 40 mK. 

We first show SB characteristics of our DQD at /3|| = 0. 
Figure 1 (b) shows a spectrum of /dot measured around 
the (l,l)-(2,0) charge boundary as a function of VpL 
and VpR. /dot is clearly suppressed inside the trapezoid 
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FIG. 1: (color online) (a) SEM image of the DQD sample, (b) 
/dot spectrum as a function of Vpr and Vpl around the SB 
regime at By = and Vc = —1.13 V. (c) Energies of relevant 
states around the SB regime at By — and t — 0. (d) Energy 
levels around e = as shown by red-boxed region in (c) for 
Bll / and f / 0. 



region at the bottom of the overlapping triangular re- 
gions (bias triangles), in which electronic transport is 
supposed to be free of CB. This suppression of /dot is 

I3I 14 [■ which originates from 



the manifestation of SB 
the energy and spin conservation relevant to the two- 
electron states shown in Figs. 1(c) and 1(d) as a function 
of the energy detuning e between the two dots measured 
from the (2,0)-(l,l) charge degeneracy point [white ar- 
row in Fig. 1(b)]. For the (2, 0) charge configuration, the 
Pauli exclusion principle makes the triplet [(2, 0)T] much 
higher in energy than the singlet [(2, 0)5'] with an energy 
gap AsT ~ 350 /ieV, making (2,0)T energetically inac- 
cessible from all other states for e < Ast- For the (1,1) 
charge configuration, both singlet [(1,1)5] and triplets 
(Ti) are within the transport window [i (— 0, ±) denotes 
the z-component (along /3||) of the total spin]. Since the 
Ti -> (2,0)5 and Ti -J> (1,1)5 transitions are prohib- 
ited by spin conservation, the sequential tunnel process 
(1, 1) ^ (2, 0) ^ (1, 0) ^ (1, 1) O for VsD < as used 
here is blocked once one of the T^'s is occupied. 

The weak features seen in the SB regime, in turn, indi- 
cate leak current that occurs only via T — )■ 5 spin tran- 
sition. In GaAs DQDs, possible mechanisms for T — )■ 5 
transition include NS-mediated processes ll^ and 
"NS-free" processes such as cotunneling [2]J and spin- 
orbit coupling [23, [2^. Contributions of these mecha- 
nisms can be distin guis hed by mapping /dot as a func- 



tion of e and /?|| 14, 3, 17 1. The data shown in Fig. 2(a) 
were obtained by sweeping e at each fixed value of /3|| . 
Consistent with the previous report for a GaAs lateral 
DQD [13], leak current appears only at e ~ (along 
the Bh axis) or at Bu ss (along the e axis, but only at 
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FIG. 2: (color online) e and /Jy dependence of /dot around 
the SB regime at Vc = —1.13 V. (a) is obtained by e sweeps 
at each value of /Jy . (b) and (c) are obtained by up and down 
sweeps of /3|| at each value of e, respectively, (d) /dot obtained 
by round sweeps of /3 1 1 around the SB regime at 1/c ~ —1.13V. 
Thick (thin) curves show the data for up (down) sweeps of 
/?y. Curves are offset vertically by 0.8 pA for clarity, (e) 
/dot obtained by very slow e sweeps for several values of /Jy . 
Curves are offset vertically by 0.2 pA for clarity. 



e > 0). As detailed in Refs. [14;l2Q[, both of these features 
originate from NS-mediated processes that take place at 
T-S crossings; the former is presumably related to the 
T_-S{j and T^-S^ double crossings that occur only at 
e « 0, while the latter involves Su and all T^'s {i — 0, ±), 
which become nearly degenerate at /3y « and e > 
[see Figs. 1(c) and 1(d)]. Here, 5y(i) denotes the upper 
(lower) branch of singlet eigenstates for finite t hybridiz- 
ing (2,0)5 and (1,1)5. Except for these regions with 
NS-mediated transport, /dot is small (less than 150fA) 
at least at e < 100 /ieV, indicating that contributions of 
NS-free processes are minor here. 

We find that strikingly different behavior shows up 
when /3||, instead of e, is swept at each fixed value of 
e. When /?y is swept up [Fig. 2(b)], a pronounced feature 
manifested as an enhancement of /dot emerges at e < 
over a wide range of e and /?j|. It is noteworthy that 
the relevant region e < is outside the bias triangles, 
which indicates that this new feature is emerging in the 
CB regime. Another important characteristic of this fea- 
ture is its strongly hysteretic behavior with respect to 
the /3y sweep direction. As shown in Fig. 2(c), the area 
of enhanced /dot becomes much smaller when /3y is swept 
down. Figure 2(d) compares /dot traces for up and down 
sweeps at several values of e. For up sweeps, the feature 
appears as a step-like enhancement of /dot , which is then 
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dragged over a wide magnetic-field range of 100 mT or 
more. For down sweeps, such dragging is much weaker, 
or even the current enhancement itself can be absent (see 
the trace for e = — 22 ^eV). 

We emphasize that the behavior of the hysteretic trans- 
port we observe at e < is distinct from that at e w 0. 
The trace for e (= —4 /xeV) in Fig. 2(d) shows that 
drag is observed for both up and down sweeps. As re- 
ported _for various material systems including not only 
^^C-enriched carbon nanotubes 
and is pre- 
sumably specific to the double crossings of T+Sl and 
T^-Su, which occur only at e sa 0. 

The observed hysteretic behavior suggests DNP in- 
duced by electron- nuclear spin flip flop ^Ifi] . The drag of 
/dot enhancement indicates that, as a result of DNP, the 
system is locked into a state that triggers the onset of /dot 
while /3|| is swept. Such locking is possible only when the 
average Overhauser field Bf^ = {Bf^ + B^ ^)/2 grows 
at a rate sufficient to cancel the variation of B\\ , a con- 
dition given by |B||| < |-B^| and sgn(B||) = -sgn(B^), 
where the overdot denotes the time derivative. Hence, 
/3^ < must hold for up sweeps {B\\ > 0), which im- 
plies DNP into spin-up because of the negative g factor 
{g* = —0.44) and a positive HF coupling constant in 
GaAs [ij]. The fact that the relevant feature appears at 
e < (and at small ) suggests the T+-Sl crossing as 
the cause of the DNP [Fig. 1(d)]. An electron spin flip 
from up to down at the T+ — Sl transition would flop 
a NS from down to up, consistent with the required sign 
of DNP. Indeed, for up sweeps the position of the onset 
of /dot behaves similarly to the calculated position of the 
T^-Sl crossing as a function of e and [dashed line in 
Fig. 2(b)], given by 2\g*\iiB{B\\ - Bo) - VSt^ + e^ + e, 
where /is is the Bohr magneton and Bq is an adjustable 
parameter. 

As shown in Fig. 2(a), DNP is absent for e sweeps at 
typical sweep rates, irrespective of sweep direction. This 
is in part due to the energy resolution of s sweeps, which 
is generally much coarser than the energy scale relevant 
to T-S mixing set by the statistical fluctuations of the 
Overhauser field ||A/3Ar[[ (typically ~ 5mT in a GaAs 



lateral DQD [14|). Indeed, for very slow and fine sweeps, 
we observe a similar /dot enhancement dragged by down 
sweeps, consistent with DNP into NS up [Fig. 2(e)]. How- 
ever, the DNP signal appears only in a narrow region and 
with much lower current level than for /3|| sweeps. We 
believe that this is because changes in the electronic state 
induced by gate sweeps are not completely canceled by 
the Overhauser field, unlike in the case of /3|| sweeps. 

We note that the T+-Sl crossing has previously been 
identified in pulsed gate experiments combined with 
charge detection, where its position as a function of e 
and /3|| has been mapped out P, @-[§]. Furthermore, 
passage through the T^-Sl crossing in a pulsed gate se- 
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FIG. 3: (color online) (a)-(d) /dot obtained by up sweeps 
of By at each value of e: (a) Vc = —1.15, (b) —1.17, (c) 
— 1.21, and (d) -1.25 V. (e) /dot obtained by e sweeps at each 
value of _B|| for Vc ~ —1.25 V. (f) Thin-solid curves show e 
dependence of c{t, e) — \ ((2, 0)S\Sl) [^ for several values of t. 
Thin-broken curves show e dependence of /?|| at the T+-Sl 
crossing point for each t. Thick-solid curves show the regions 
with c{t, e) > 10'^ at e < 0. 



quence has been harnessed to dynamically polarize nu- 
clei [1, 0] or to prepare a preferred NS state @. In 
DC transport, on the other hand, due to CB at £ < 0, 
which prohibits the charge transfer from (1,1) to (2,0) 
involved in the sequential tunnel process, the T^-Sl 
crossing is not usually observed [25]. We suggest that 
in the present case CB is lifted by making the tran- 
sition from (1,1) to (1,0) through Sl, that is, with- 
out making a real transition to (2, 0)5". The transition 
from Sl to (1,0) is made possible by a finite (2,0)5* 
component in Sl that exists even at e < for finite 
t. The amplitude of this (2, 0)5" component is give n by 
cit,s) ^ mo)S\SL)\' = [1 + (eA)/yW+8]/2 H, 
indicating that it decreases with increasing \e\ at e < 0. 
This qualitatively explains the observation that the DNP 
signal decays with increasing \e\ [Fig. 2(b)]. 

The above equation shows that c{t, e) is a function of 
e/t, suggesting that the ratio e/t, rather than e itself, 
is the key parameter dictating the occurrence of DNP. 
Figures 3(a)-3(d) show /dot obtained by up sweeps of /?|| 
for several Vc and thus for different values of t. While 
a DNP signal is present for Vc > —1.21V, as Vc is 
made progressively more negative, the area of e where 
DNP is observed shrinks. By Vc = -1.25 V, DNP is 
no longer discernible, where the overall feature becomes 
nearly identical to that for the e sweeps [Fig. 3(e)]. The 
thin solid lines in Fig. 3(f) depict c{t,e) as a function of 
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FIG. 4: (color online) Compositions of (2, 0)5, (i, t), and (f, i 
) states in (a) S'j^ and (b) Tq as functions of \g*\fiB^B^. e — 
—40 ^eV and t = 3 /ieV are assumed. Grey regions indicate 
the range of statistical fluctuations, jAB^I < ||Ai3jv||/2. 



e for several values of t. Since smaller t gives smaller 
c{t,e) at e < 0, the range of e with a certain value of 
c{t,e) shrinks with decreasing t. Figure 3(f) also shows 
the Tj^-Sl crossing point in the plane of e and Bn for [ll|, [l 
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several values of where regions with c(t, e) > 10 
are highlighted with thick lines. These regions become 
shorter with decreasing in qualitative agreement with 
the observed e and B\\ dependence of the onset of DNP 
[Figs. 2(b) and 3(a)-3(c)]. 

It is noteworthy that, different from the previous re- 
port on a vertical DQD [ll], DNP due to the T_-Su 
crossing at e > is not observed here. Note that DNP via 
the T--Su crossing is predominantly limited by the tran- 
sition rate from T+ to other states, which are higher in 
energy at a finite . Since this requires phonon absorp- 
tion for energy conservation, it is strongly suppressed at 
low temperature in the absence of NS-free processes such 
as cotunneling, as is the case here. This contrasts with 
the case of the T^-Sl crossing, where DNP is mainly lim- 
ited by the relaxation rate from T_ to lower lying states 
via HF coupling, which can take place with phonon emis- 
sion even at low temperature. 

While the occurrence of DNP and its dependence on 
£ and t can be explained as above, the dependence of 
/dot on B|| provides a deeper insight into the underly- 
ing transport mechanism. That is, when the dot state 
is locked via DNP, /dot does not remain constant, but 
increases with [Fig. 2(d)], which points to the exis- 
tence of an additional degree of freedom that governs 
the transport. We propose a model that DNP generates 
imbalance ABf^ in the NS polarization between the two 
dots and this finite ABf^ enhances /dot- When ABf^ ^ 0, 
Shi Su, and Tq are no longer eigenstates of the system, 
and mix to form new eigenstates 5'^, S'jj, and Tq. Fig- 
ures 4(a) and 4(b) depict how the characters of S'j^ and 
evolve with AEz = \g*UiB^BI, using (t,!), a,t), 
and (2,0)5' as a basis set ^q. Here, (t,|) [(|,t)] de- 
notes the state with spin-up (-down) and spin-down (- 
up) electrons in the left and right QDs, respectively. The 



calculation shows that, for AEz > 0, the (t,i) compo- 
nent in S'j^ becomes progressively dominant with increase 
in AEz- This implies that for AEz > the probabil- 
ity of S'j^ transition being mediated by a NS flip 
in the right QD increases with AEz- The key observa- 
tion here is that the resultant down-to-up NS flip in the 
right QD increases AEz and further facilitates NS flip in 
the right QD. The same argument holds for AEz < 0, 
which facilitates NS flips in the left QD. Due to this pos- 
itive feedback, large |A/?|_f| can grow spontaneously out 
of tiny statistical fluctuations. This flnite ABf^ produces 
a dramatic consequence on transport: when AB^ ^ 0, 
Tq has a finite (2, Q)S component, which lifts the SB due 
to Tq. With increase in |A/3^|, this (2,0)5" component 
first grows rapidly at small |A/3^|, but then increases 
only gradually at large |A/3^|, which accounts for the 
observed sharp onset and subsequent gradual increase of 
/dot • Imbalance in NS polarization has also been invoked 
to account for the hysteresis and multiple-step behavior 
of leak current observed in a DQD under ESR condition 
Despite the marked differences in the experi- 



mental situations between the present study and that in 
Ref. [lH, AB^ plays critical roles in both experiments, 
which clearly demonstrates that it must be taken into ac- 
count to fully understand electron-nuclear spin dynamics 
in DQDs, not only in pulsed-gate, but also in DC-drive 
experiments. 

In summary, we have reported a new hysteretic trans- 
port feature in the CB regime of a GaAs lateral DQD, 
which signals DNP due to the T+-Sl crossing. This DNP 
signal appears owing to (2, 0)5-hybridization in Sl by fi- 
nite t- The magnetic-field dependence of the current level 
in the DNP region suggests strongly unbalanced NS po- 
larization. 

We thank Y. Okazaki and T. Ota for experimental ad- 
vice and Y. Tokura and Y. V. Nazarov for discussions. 
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HAMILTONIAN 

Eigen energies are calculated with the Hamiltonian H = -ffoQD + H^, where -ffoQD and are 
Fdqd = -e\{2,0)S}{{2,0)S\+ J] ag>BB|||r,)(T,| 

cr = ±l 
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on a basis set of (2,0)5, (t,i), (i,t), T+, and T_. Note that |(t,i)) = [\Tq) + 1)5)]/V2 and |(i,t)) = [|To) - 
|(1,1)5)]/V2 for derivation of Figs. 1(c), 1(d) and 3(f). Fi gurelST] shows typical energy level diagram for the five eigen 
states as a function of g* fiB^Bf^. Since only the 5^ level can cross at e < with the T+ level from this figure, only 
the r+ — > S'j^ scattering is considered in our letter. 
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FIG. SI: Energy level diagram of the five five eigenstates 5^, Tq, S^, T+, T— for the H as a, function of g* fiB^B^ . e ■ 
and g*fiB{B^\ + B^) = —LOfieV are assumed. 
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